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ABSTRACT 

Context. Recent high-resolution observations have shown that stellar winds harbour complexities that strongly deviate from spherical 
symmetry, which generally is assumed as standard wind model. One such morphology is the Archimedean spiral, which is generally 
believed to be formed by binary interactions, as has been directly observed in multiple sources. 

Aims. We seek to investigate the manifestation in the observables of spiral structures embedded in the spherical outflows of cool stars. 
We aim to provide an intuitive bedrock with which upcoming ALMA data can be compared and interpreted. 

Methods. By means of an extended parameter study, we modelled rotational CO emission from the stellar outflow of asymptotic 
giant branch stars. To this end, we developed a simplifled analytical parametrised description of a 3D spiral structure. This model 
is embedded into a spherical wind and fed into the 3D radiative transfer code LIME, which produces 3D intensity maps throughout 
velocity space. Subsequently, we investigated the spectral signature of rotational transitions of CO in the models, as well as the 
spatial aspect of this emission by means of wide-slit position-velocity (PV) diagrams. Additionally, we quantifled the potential for 
misinterpreting the 3D data in a ID context. Finally, we simulated ALMA observations to explore the effect of interefrometric noise 
and artefacts on the emission signatures. 

Results. The spectral signatures of the CO rotational transition v=0 J=3-2 are very efficient at concealing the dual nature of the 
outflow. Only a select few parameter combinations allow for the spectral lines to disclose the presence of the spiral structure. If the 
spiral cannot be distinguished from the spherical signal, this might result in an incorrect interpretation in a ID context. Consequently, 
erroneous mass-loss rates would be calculated. The magnitude of these errors is mainly confined to a factor of a few, but in extreme 
cases can exceed an order of magnitude. CO transitions of different rotationally excited levels show a characteristical evolution in 
their line shape that can be brought about by an embedded spiral structre. However, if spatial information on the source is also 
available, the use of wide-slit PV diagrams systematically expose the embedded spiral. The PV diagrams also readily provide most 
of the geometrical and physical properties of the spiral-harbouring wind. Simulations of ALMA observations prove that the choice 
of antenna configuration is strongly dependent on the geometrical properties of the spiral. We conclude that exploratory endeavours 
should observe the object of interest with a range of different maximum-baseline configurations. 

Key words. Line: proflles-Radiative transfer-Stars: AGB and post-AGB-circumstellar matter-Stars: winds, outflows- 
Submillimeter: stars 


1. Introduction 

Near the end of their lives, low- and intermediate-mass stars as¬ 
cend the asymptotic giant branch (AGB), where they develop a 
dense molecular envelope interspersed with (sub-)microscopic 
dust particles. The rate at which these stars expel their outer 
envelopes is high, from a few times 10 ~^ M© yr“^ to about 
10“"^ M© yr“^ (e.g. iDe Beck et al.l 1201 d), and often va ries in 
strength on a range of timescales (e.g. lOlofsson et al.l I201QI: 
iMaercker et al.l [201^ . The winds are thought to be driven by 
a combination of pulsations and radiation pressure on dust 
grains, although the fundamentals of this mechanism have 
turned out to be an elusive problem. So far, empirical mod¬ 
els and hydrodynamic simulations of these outflows have fo¬ 


cused mainly on describing isolated AGB stars featuring spher¬ 
ical outflows. Recent observations with high-resolution tele¬ 
scopes, however, have revealed a rich spectrum of structural 
comp lexities. T hese include bipolar structures (e.g.fealick et ^ 
20131). arcs te.g.lPecin et al.ll2012l:ICox et al.ll2012h. shells fe.g. 


Mauron & Huggins |200Q|). clumps (e.g. iBowers & Johnstod 
199Qh. spirals (e.g.lMauron & HugginJ2006l:lMaver et al.l201ll: 

ar " 
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Maercker et al.ll2Q12l: iKim et al.ll2013l). tor i ife.g. ISkinner et alJ 

19981) . and bubbles (iRamstedt et al.l 120141) . Several processes 
may be responsible for the formation of these structures. Shells 
and arcs may be caused by tem poral variations in the mass-loss 
rate and/or expansion velocity (IMaercker et al .11201^ and refer¬ 
ences therein). Aspheri cal structure s may be the result of non¬ 
isotropic mass-loss fe.g. lUetall2006), systemic motion relative to 
the local interstellar medium (^e.g. iDecin et aLllMl^ . magnetic 
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fields fe. g.lPerez-Sanchez et al.ll20131: Ivan Marie et al.l[2014l) . or 
binarity dSokerll 1 9971: iHugginsll2007h . The latter need not be sur¬ 
prising as the multiplicity frequency of the progenitors of AGB 
stars is above 50 percent (iRaghavan et aLlSoTO : iDuchene et al.l 
l2Q13l) . A complex morphology of the circumstellar envelope 
(CSE) may have an effect on the strength and shape of spec¬ 
tral lines as well. Exotic profile shapes may result from devia¬ 
tions from s phericity, which may not be recovered using spher¬ 
ical models (iDe Beck et al .11201^ . Alternatively, the lines may 
have shapes typical for spherical flow but strengths that are de¬ 
viant. In the latter case, spherical models may lead to erroneous 
determinations of the overall wind properties. 

In the present paper we attempt to broaden our under¬ 
standing of embedded stellar wind morphologies by develop¬ 
ing simplified mathematical models of stellar wind structres. 
Here, we focus on th e spiral structure, for which very convinc 


2012 

), and theoretical background exists (iTheuns & Jorissen 

1993 

: ISoked 11994 Mastrodemos & Morris! 119991: 

Kim 2011: 

Mohamed & Podsiadlowski 2012). Hvdrodvnamical Simula- 


tions indeed show that the influence of a binary companion on 
its direct surroundings can form Archimedean spirals embedded 
in the outflow, along the orbital plane of the system. A few dis¬ 
tinct mechanisms are known to produce such a morphology. The 
existence of a binary companion inside the CSE of the mass¬ 
losing star can cause local perturbations which, after propaga¬ 
tion through the CSE, develop spiral patterns. Local companion- 
wind interactions cause the wind material near the companion 
to be gravitationally funneled into an enhanced density region (a 
gravity wake) resulting in a flattened (i.e. only barely extending 
away from the orbital plane) spiral. Additionally, the presence of 
this binary companion will cause the mass-losing star to wob¬ 
ble around the common centre of mass of the binary system. 
This reflex motion alters the local densities, ultimately forming 
an Archimedean spiral that is comprised of conce ntric shells , 
when observed edge-on (for visualis ation, see EigO (lKimll2Qlll: 
iKim & Taamll201^lKim et al]l2013l) . Here, the specific geomet¬ 
rical properties of the resulting spiral have been found to be de¬ 
termined by the strengt h of the local pavit y field with respect 
to the outflow velocity dKim & Taamll2012b . Radial models of 
shocks generated by a pulsating mass-losing source, that are lo¬ 
cally enhanced by interaction with a binary companion have also 
shown that not on ly single, but also multi ple spiral branches can 
form in the CSE (IWang & Willson|[2012h . In this case the prop¬ 
erties of the spiral depend strongly on the ratio between the pul¬ 
sation period of the star and the binary period. Einally, if the 
binary contains two mass-losing stars, the c olliding winds at the 
interaction zone may form a spiral pattern (IStevens et al.lIl992t 

IWalder & Eolini|[200^ Ivan Marie et al. 11201 ik . Here the geomet- 

rical properties of the spiral are determined primarily by the ve¬ 
locity of the slower wind and the particular orbital motion of the 
binary. 

In this paper, we make no assumptions on the spiral forma¬ 
tion mechanism. Instead, we assume its presence, and show how 
these structures are manifested in the observables by means of 
3D radiative transfer calculations. The stellar wind properties are 
described analytically, and its emission is modelled with LIME, a 
fully three-dimensional non-local -thermodynamical-equilibrium 
(NLTE) radiative transfer code (iBrinch & Hogerheijdel 120101) . 
Arguably, an exclusively analytical description of the stellar 
wind properties might not be very physically consistent. How¬ 
ever, this approach allows us to calculate a considerable grid of 
models, which at the present time would be virtually impossible 
to achieve via consistent radiative-hydrodynamical calculations. 


We present the results both in spectral form, exhibiting the man¬ 
ifestation of the spiral structure in the molecular emission lines, 
and in spatial form, providing a reference for high spatial resul- 
tion data from telescopes like ALMA. 

This paper is organised as follows: In Sect. 2 we describe the 
mathematical structure of the used analytical expression in de¬ 
tail, followed by the computational tools which have been used 
to produce the results. In Sect. 3 we give an overview of the ex¬ 
plored parameter space, and of the assumptions we have made 
in this paper. Section 4 presents the results of the radiative trans¬ 
fer, both in spectral and in spatial format. Einally, in Sect. 5, 
we make an effort to relate the emission to intrinsic geometrical 
properties, and we show the effect of ALMA antenna configura¬ 
tions on the characteristics of the intrinsic emission. 


2. Numerical procedure 

2.1. Spiral geometrical model 

The following mathematical descriptions are formulated in the 
spherical coordinate system, 

r = ^x^+y^+z^ (1) 

6 = arccos^-j (2) 

0 = arctan^-j, (3) 

representing radial, azimuthal and equatorial distance respec¬ 
tively. The Archimedean spiral is the solution of the equation 
r = b(f), and is characterised by a constant distance (d = 2nb) 
between consecutive turnings. The analytical properties of the 
assumed spiral model are described by a Gaussian distribution, 
such that sharp cut-offs of the wind properties described by this 
geometrical pattern are avoided. In dimensionless units this dis¬ 
tribution is given by 


S (r, 6, (p) = exp 


(r - ro(r, 6>, (p)f _ {0 - 9o(r, 0, 
2(Tr{r, 0, (pf 2cr(,(r, 0, (f>)^ 


(4) 


where and ctq are Gaussian widths in radial and azimuthal 
space. In the following subsections, we describe the specific ge¬ 
ometrical properties of this model both radially and azimuthally, 
and provide an in-depth explanation of the functional parame¬ 
ters present in Eq. IH The connection between S (r, 6, 0) and the 
distributions of stellar wind properties (density, temperature) are 
outlined in Sect. 3. 


2.1.1. Geometry in the orbital plane 

The parameters describing the geometry in the orbital plane are 
given below and are illustrated in Eig.[Tl 

- ro(r,0,(p) = bicp - 0o) describes the exclusive 0-dependent 
outward evolution of the spiral windings. With 0 > 0 and 
00 > 0, 00 corresponds to the position angle, that is, the 
angle with respect to a reference point from which the spiral 
originates. 

- cTrir, 6, 0) = cTr describes the Gaussian width of the spiral 
arm in the orbital plane. 
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Fig. 1: Face-on view of the spiral model by taking a slice through 
the orbital plane, qualitatively depicting the meaning of the geo¬ 
metrical parameters. The colour scheme used is qualitative, and 
serves to demonstrate the Gaussian nature of the S (r, 6, (p) func¬ 
tion. 

2.1.2. Geometry in the meridional plane 

The parameters describing the geometry in the meridional plane, 
which is perpendicular to the orbital plane, are: 

- 6Q{r, 6, (p) = 7r/2, chosen such as to keep the Archimedean 
spiral symmetrical relative to the orbital plane. 

- (T6i(r, (p,6) = a determines the angular height of the spiral, 
perpendicular to and symmetrical around the orbital plane. 
It is, in effect, a measure for the extent by which the spiral 
fans out away from the orbital plane. Its mathematical equiv¬ 
alence to the Gaussian standard deviation in Eq. |4] implies 
that, in case of a density spiral, 68.2% of the total mass in 
the spiral can be found within the volume bound by a. The 
choice of an angle as the parameter to quantify the height 
of the spiral is justified by the geometrical implications of 
an exclusively radial outflow. If a is small, then the bulk 
of the material is confined into a narrow spiral close to the 
orbital plane. If caused by binary interactions, this would 
correspond to a case where the spiral is formed by a strong 
gravitational wake of the companion. For larger a, the spiral 
windings become more extended shell-like structures. Such 
a morphology would be a reasonable representation of the 
effect of a wobbling mass-losing star. In the most extreme 
case, which is for an opening angle of a = 180°, one is left 
with a dense winding of shells. A visualisation of the effect 
of a on the edge-on view of the spiral in shown in Fig.[^ 


2.1.3. Additional considerations 

Hydrodynamical simulations bv iKimI (1201 ll) show the existence 
of morphological substructures within the overall spiral shape. 
These models show this substructure to consist of two adjacent 
spirals. When seen edge-on (as in Fig.O a frontal (curved) wake, 
and a rear (flat) wake can be identified. Such details can also be 
modelled with the previously developed analytical expression. 
The frontal spiral can be modelled by adopting Eq. |4l To model 
the rear portion of the spiral substructure, Eq.|4]needs to be trans¬ 
formed from spherical coordinates to cylindrical coordinates, by 



Eig. 2: Cross-section through the centre of, and perpendicular 
to the orbital plane. The effect of an increase in a on the edge- 
on view of the spiral model, with a the angle between the two 
dashed lines. A narrow spiral with an opening angle of a = 10° 
(left panel), an intermediate spiral with opening angle a = 35° 
(middle panel), and a shell spiral with a = 180° (right panel) are 
shown. The red horizontal dotted line represents the location of 
the edge-on orbital plane. The colour scheme used is qualitative, 
and serves to demonstrate the Gaussian nature of the S (r, 6, (p) 
function. 

substituting the spherical radius (r = with the 

cylindrical radius (r = a/?T^), and the spherical height (6) 
with the cylindrical height (z). By subsequently correctly choos¬ 
ing the ratio between the Z^-parameters of both spiral compo¬ 
nents, a resulting spiral pattern consistent with hydrodynamical 
results can be constructed. 

We recall the probability of observing the orbital plane of a 
binary system face-on or edge-on. We define the face-on position 
to be an orientation for which the normal on the orbital plane de¬ 
viates by at most ten degrees from the line of sight. Similarly, 
we define the edge-on position of the orbital plane to be the ori¬ 
entation for which the normal deviates by at most ten degrees 
from being perpendicular to the line of sight. With these defi¬ 
nitions, we find that the probability of a face-on orientation is 
1.5%, compared to 17.3% for edge-on. One is thus more than 
ten times more likely to observe a binary system edge-on than 
face-on. 

We also make a note on our definition of position angle, 
which, in effect, is a rotation around the z-axis (or a rotation 
over the angle (p). However, the rotational operations over 6 and 
(p are not commutative. If the spiral is rotated over (p first, then 
the Archimedean spiral will remain in the x,y-plane. Its effect 
will be a face-on rotation, revolving the start of the spiral around 
the origin of the coordinate system. This effect is trivial, and is 
ignored in this paper. If, on the other hand, the spiral is rotated 
around 6 prior to a manipulation of 0, then the inclined spiral 
will be pivoted around the z-axis. Under such an operation, the 
normal on the Archimedean spiral plane will precess around the 
z-axis. We refer to angular manipulation of the latter case as op¬ 
erations on the position angle. 

2.2. Radiative transfer 

To generate the 3D intensity channel maps of simulated 
emission of the analysed spiral models, the Non-Local- 
Thermodynamica l-Equilibrium (NLTE) full-3p radiative trans¬ 
fer code LIME (iBrinch & Hogerheijdel 120101) was used. Eor 
a technical overview of the inner workings of the code, see 
iBrinch & Hogerheijdel (1201 Oh . 

The code was originally designed for modelling cool molec¬ 
ular clouds. Eor the modelling of stellar winds, however, the cen¬ 
tral mass-losing star is a non-negligible source of energy. We 
have implemented the option to position an arbitrary number of 
stars at locations of choice. Radiative interaction of every grid 
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point with these ’stars’ is forced at least once per iteration. The 
stars emit as black bodies. It is important to note that the code 
is non-dynamical, meaning it exclusively calculates the radiative 
transfer for the input parameters. 

We modelled CO emission, for which the spectr oscopic CO 
data of the LAMDA database (ISchoier et al.ll2005[) we re used. 


The collisional rates were taken from 


YangetalJd2m(l . 


2.3. Synthetic ALMA simulations 

To produce simulations of ALMA observations of the emis¬ 
sion of the spiral wind models we used the Common As- 
tronomy Software App lications, CASA, post-processing package 
jMcMullin et al.|[200% . It consists of a collection of C++ tools, 
managed by a python wrapper code. 

The specific tools used to simulate the observations are the 
‘simobserve’ and ‘clean’ functions. The simobserve task con¬ 
verts any model image (corresponding to the true sky brightness 
distribution) into virtual observations. For that it is necessary to 
specify quantities such as integration time, antenna configura¬ 
tion, frequency, and astrometrical and atmospherical parameters. 
These are used to calculate the noise contributions in uv-space. 
Subsequently, the task creates the Fourier transform of the model 
image and projects this onto this grid. The clean task converts 
your visibilities from uv-space back into real coordinates, to sim¬ 
ulate how ALMA would see the input model. Finally, it performs 
a deconvolution of this ‘dirty image’ to get rid of the sidelobe 
structure, yielding a 3D datacube of a synthetical ALMA obser¬ 
vation of the intrinsic 3D emission. 


3. Model assumptions 

We consider a singular mass-losing system, located in the centre 
of the model. Though likely to be formed by binary interactions, 
we attempt to make no assumptions on the formation of the spiral 
structure in this paper. For all intents and purposes, this singular 
mass-losing system represents the binary system. 

The 3D density field consist of two components. One spheri¬ 
cal outflow component, for which mass conservation implies that 

PhoC'*) = MHo/(47rrh(r)). (5) 

Mho is the rate at which mass is lost in the homogeneous out¬ 
flow. Hydrodynamical models show the outflow v elocity of the 
spiral to be compyable to the intrinsic wind speed (iKim & TaamI 
1201 2l: iKim et al.] l2Q13l) . We express the velocity of both the 
spherical and the spiral wind components by a beta-law, 

v(r) = Vc^(l -Ro/rf, (6) 

where Vc» is the terminal wind velocity, and Rq the dust conden¬ 
sation radius. v(r) is thus exclusively radial in nature, and is the 
key to understanding why the spiral is bound by the opening an¬ 
gle a. In addition to v(r) a constant turbulent velocity field v^urb 
exists throughout the wind. 

This homogeneous density field is superimposed by a spiral 
density enhancement. 


Pspirai(^, 0) = Po(l 6, 0) S (r, 6, 0) (7) 

where S (r, 6, 0) is described in Eq. |4l po(r, 6, 0) = !r)^ 

describes the variation of the local Gaussian density peak, where 
Pmax is a constant, and will be determined later. Supported by 
mass conservation in a radial outflow, we assumed a radial 
dependence of po proportional to 1/r^. The resulting overall 


density profile of the spiral wind is Pwind(^, 0) = Pho(^) + 
PspiralC^? ^9 0)- 

The temperature distribution throughout the wind follows the 
same twofold format as the density structure. Its main trend is a 
radial power law, 

Twir) = T.(R./rr, ( 8 ) 

which describes the temperature of the homogeneous outfiow. 
Thc and R^ are the stellar effective temperature and stellar radius 
respectively. The temperature of the spiral, 

7^spirai(^, 0, 0) = To(r, 6, 0) S (r, 6, 0), (9) 

is superposed onto this background, with a chosen distribu¬ 
tion maximum of To(r, 6, (p) = Tq = 60K, corresponding to a 
value of the same order of magnitude as estab lished from hy¬ 
drodynamical calculations (iKim & TaamI l2012h . The resulting 
overall temperature profile of the spiral wind is T^mdiL 0, 0) = 
7^Ho(^) + 7^spiral(^, 0 , 0). 

Though not expected to affect the CO emission much, dust 
was taken into account in calculating the radiative-transfer, with 
the dust density distribution following the gas density distribu¬ 
tion. The specific dust input parameters are presented in Sect. 
2.2.4. 

Additional effects such as scattering, magnetic fields, or ro¬ 
tation are not taken into account. 

3.1. Outflow Parameters 

We conducted a parameter study to asses the effect of changes in 
the free parameters of our model. Here, we present the applied 
parameter values. The general stellar and dust properties, as 
well as some overall characteristics of the CSE were taken from 
iDe Beck et al.l (l2012h . making the high mass-loss case of the pa¬ 
rameter study similar to the carbon-rich CW Leo system. The 
bias towards C-rich objects reflects that most spirals were ob¬ 
served in such environm ents (iMorris et al.ll2006l:lMaercker et al.l 
I2ni2l:lnecin ^1^120151) . 


Stellar Parameters 


n 

2330 K 

L, 

11300 Lo 

Distance 

150 pc 

Wind Parameters 

Voo 

14.5 km s ^ 

^turb 

1.5 kms“^ 

P 

0.4 

Ro 

2.0 

€ 

0.5 

CO/H2 

6.0 X 10-4 

Spiral Parameters 

Inb 

270 AU 

(Tr 

20 AU 

To 

60 K 

Dust Parameters 

Amorphous Carbon 

53% 

Silicon Carbide 

25% 

Magnesium Sulfide 

22% 

Gas/Dust 

100 


Table 1: Fixed model parmeters. 

Table [D gives an overview of the fixed parameters of the 
radiative-transfer models. These are the parameters which are 
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either deemed to have a predictable or even trivial effect on the 
observables, or which bring about a global influence which is 
beyond the scope of this paper. The fixed spiral parameters are 
chosen such that a spiral with sufficiently contrasting features 
can be clearly identified in the observables. 

The variable parameters are presented in Table [H and the 
motivation behind their choice is described below. 

- The total mass-loss of the system has profound effects on 
the emission in terms of absolute strength and optical depth 
effects. Therefore we simulated CSEs for two extreme mass- 
loss rates, one very high and one very low. 

- The contrast between the density of the spiral and the homo¬ 
geneous outflow brings about changes in relative emission 
strengths, which can easily be recognised in the observables. 
We define a mass contrast between both wind elements, us¬ 
ing the parameter 


oriented such as to maximally exploit the asymmetry of the data. 
This also counteract additional projection effects brought about 
by the position angle of the system. Our PV diagrams were con¬ 
structed following these instructions. The axes along which the 
asymmetry of the models appears strongest are labelled as X and 
Y. Much easier to interpret than the full 3D data, they provide 
the user with clearly correlated structural trends of the complete 
wind, which in turn provides strong clues on possible geometries 
harboured by the wind. 

All the shown data will are of the emission of the ground vi¬ 
brational CO rotational transition J=3-2, unless otherwise stated. 
All results were continuum subtracted. In this section, the PV di¬ 
agrams that best exhibit the asymmetry of the data (along X and 
Y) are referred to as PV 1 and PV2 respectively. The spectral data 
were modelled as seen with a 20 arcsecond beam. (The quality 
of most images shown below is strongly enhanced when viewed 
on screen.) 


^ _ ^spiral _ ^spiral (10) 

^HO ^HO 

with Mspirai + ^HO = Motai- ^spiral represents the portion 
of the material confined to the boundaries of the spiral struc¬ 
ture, flowing outward per unit time. S thus quantifies the way 
in which the total mass lost by the star is distributed over 
the spherical outflow and spiral density enhancement com¬ 
ponents. S ^ 0 corresponds to Mho ^ MotaU S ^ oo to 
^spiral ^ Motai- A particular choice of S will result in a spe¬ 
cific value for Pmax- In principle, this parameter should also 
affect the value of Tq. We have chosen not to let it depend on 
E because the relation between both is not clear to us. 

- Two different spiral heights were simulated: one narrow spi¬ 
ral, with a small opening angle, and one spiral having a max¬ 
imally large opening angle, referred to as a shell spiral. 

- The inclination i is defined as the angle that tilts the view of 
the model between face-on, where the orbital plane lies per¬ 
pendicular to the line of sight, and edge-on, where the orbital 
plane lies along the line of sight. Each model is observed un¬ 
der six different inclination angles, evenly spaced between 
this face-on and edge-on view. 


4. Results 

We present the general findings of the parameter study in two 
parts. Eirst, we exhibit the general trends and effects of the pa¬ 
rameters on the spectral lines. Second, we show the effect of 
identical parameter simulations on the intensity maps. In order 
to investigate these 3D data we make use of the so-called posi¬ 
tion velocity (PV) diagram. The three dimensional data consists 
of two linearly independent angular dimensions, representing the 
two angular coordinates in the plane of the sky, and one veloc¬ 
ity dimension. The PV diagram is, in effect, a slice through the 
3D data at an arbitrary angular axis in the plane of the sky (thus 
preserving the velocity axis). It is thus a 2D plot of the emis¬ 
sion along this chosen axis, versus velocity. In principle any slit 
width can be chosen. If the slit width is larger than one singu¬ 
lar data pixel, then the emission is collapsed onto each other by 
summing up the emission with identical PV coordinates. Eor a 
condensed overview of the overall wind morphology, it is useful 
to make PV diagrams with a maximal slit width, namely the full 
size of the datacube. This is referred to as a wide-slit PV dia¬ 
gram. Additionally, it is advantageous to construct two different 
wide-slit PV diagrams of the 3D data, choosing any set of lin¬ 
early independent (and thus perpendicular) angular dimensions 


4.1. Reference model 

Pig. [3] shows both the spectral line and the PV diagrams of 
an optically thin narrow spiral, with a total mass-loss = 1.5 x 
10“^ Mo yr“^ and 'Z = 1. This is considered the reference 
model, and will be referred to as such. Every model below shows 
the effect of the change of a singular parameter with respect to 
this model. 





0 max 


Pig. 3: Top Left: The spectral line of the reference model (blue), 
seen from an inclination angle of 0 degrees (face-on). Overplot¬ 
ted is the spectral line of an exclusively homogeneous spheri¬ 
cal outflow with identical global characteristics. Top Right: The 
wide-slit PV diagram of the reference model. The linear colour 
map at the bottom provides a reference to the colour coding. Ev¬ 
ery coloured diagram is constructed using this map. MAX rep¬ 
resents the maximum emission in the PV diargam. 


4.2. Spectral aspect 
4.2.1. Narrow spiral 

Inclination: Pig. 0] shows the effect of the inclination angle un¬ 
der which the spiral is seen. The appearance of the narrow spiral 
is very dependent on its orientation with respect to the observer. 
We therefore expect the inclination to have a strong effect on the 
observables. The homogeneous background wind is indifferent 
to these same transformations. The face-on view (i = 0) of the 
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Parameter 

Values 

Fabels 

Mass loss M [Mo/yr] 

1.5x 10-'> 

1.5x 10-’ 

H (‘High’) 

L (‘Low’) 

Mass contrast Z = nim 

V ^e{l}:m€{2,5,10,100} 

V me {1} : ^6 {1,2,5,10,100} 

Sd2,Sd5,Sdl0,Sdl00 

S1,S2,S5,S10,S100 

Spiral height a 

o 

o O 

O 00 

N (‘Narrow’) 

S (‘Shell’) 

Inclination 6 = (^/5)(7r/2) 

^6 10,1,2,3,4,5} 

1=0,18,36,54,72,90 


Table 2: Variable model parameters. The labels are used throughout the paper to specify the considered models. 


special line clearly shows the dual nature of the CSE. The spi¬ 
ral wind produces a sharp and narrow peak in velocity space. 
The width of this peak depends strongly on the height of the spi¬ 
ral. It is much narrower than the terminal velocity of the wind. 
This is because the spiral only hardly extends away from the 
V = 0 km s“^ Doppler plane. Because the wind is optically thin, 
higher local densities generate more emission, which is why the 
peak extends above the lower plateau, generated by the spherical 
wind. As the inclination angle increases, the contribution of the 
spiral decreases and widens, whilst the contribution of the ho¬ 
mogeneous background wind remains unchanged, as expected. 
This evolution of the central peak is explained by the fact that 
the total emission is smeared out over a wider velocity range. At 
the highest inclinations, a double-peaked profile emerges in the 
spiral wind feature. Additionally, its width becomes comparable 
to the width of the spherical outfiow, ultimately concealing the 
dual nature of the system. As a side note, the position angle has 
no effect on the spectral lines. 

Mass contrast: The main effect of the mass contrast is on the 
relative heigths of the spectral features generated by the spiral 
and spherical winds, as seen in Fig. [5] For E > 1, the contribution 
of the spherical wind rapidly diminishes and disappears from the 
resultant profile, and only the spectral feature generated by the 
spiral wind remains. For E ~ 1 we find that the relative sizes of 
both contributions are comparable. When E becomes very small, 
the spherical wind dominates the spectal feature. 


i = 0 i = 90 



Fig. 5: Spectral lines of a narrow spiral in the low mass-loss 
wind, in which the total amount of mass contained in the homo¬ 
geneous outfiow is five times the amount in the spiral (Z = 1/5), 
seen face-on (left panel) and edge-on (right panel). 

Figure [5] shows a model with Z = 1/5. The face-on case 
shows both a decreased contribution of the spiral wind and an 
increased contribution of the spherical wind relative to the ref¬ 
erence model. Viewed edge-on, the spiral wind’s double-peak 
characteristics (visible in theZ = I case) are reduced, leaving a 
fiat-topped resultant spectral profile. 

Total mass-loss: The first and most important effect of an in¬ 
crease in the total mass-loss is an increase in flux, visible in Fig. 
[6j For low-inclination angles the peak due to the contribution 
of the spiral wind broadens with respect to the face-on reference 


model. This is due to the increased emission in the Gaussian tails 
of the density distribution of the narrow spiral. For high inclina¬ 
tions (nearing an edge-on view) the double-peaked aspect of the 
spectral line is less pronounced as the spiral becomes optically 
thick and relative brightening due to its long column depths in 
the line of sight direction (at maximum and minimum velocity) 
is suppressed. 


i = 0 i = 90 



Fig. 6: Spectral lines of a narrow spiral in the high mass- 
loss wind, in which the total amount of outstreaming mass is 
equally divided between both wind components, seen face-on 
(left panel) and edge-on (right panel). 


4.2.2. Shell spiral 

In the case of the shell spiral, the spiral and spherical wind 
components possess virtually identical widths in velocity space, 
which makes both components essentially indistinguishable 
from one-another. This is due to the broadening effect of the 
opening angle a on the emission feature of the spiral wind. For 
high values of a both wind components will be of comparable 
width. For this reason we refrain from showing this in figures. In 
Sect. 4.3 we discuss the integrated line fluxes for various com¬ 
binations of mass-loss in spiral and spherical outfiow to assess 
the potential for erroneous estimates of the total mass-loss rate 
if assuming a spherical wind only. 

4.2.3. Other CO lines 

Different rotational transitions of CO (besides CO J=3-2) can 
provide additional information. In Fig. [T] the effect of rotational 
CO transitions ranging from J=3-2 to 1=38-37 (probing the in¬ 
ner parts of the wind) on the line shapes of the reference model is 
shown. Hydrodynamical simulations of the regions close to the 
mass-losing star show that the sphericity of the outfiow is com¬ 
pletely destroyed locally by the binary interactions. The effect 
of a well-behaved analytical spherical outfiow close to the cen¬ 
tre of the spiral is thus unrealistic. The spiral shape, however, is 
recognisable on this scale. Therefore, to properly show how the 
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i = 0 i = 18 i = 36 i = 54 i = 72 i = 90 



Fig. 4: Spectral lines of a narrow spiral in a low mass-loss wind, in which the total amount of outstreaming mass is equally divided 
between both wind components, as seen from a range of different inclination angles. 


spectral lines are affected by the spiral, we omitted the contribu¬ 
tion of the spherical wind. 


J=3-2 J = 10-9 J = 17-16 J=24-23 J=31-30 J=38-37 



Fig. 7: The effect on the shape of the spectral lines as a function 
of CO transition. The first row shows the line shapes for a face- 
on narrow spiral, the second row for an edge-on narrow spiral, 
and the third row for a shell spiral. To focus the attention on the 
shapes of the emission lines, the peak strengths have been nor¬ 
malised, and the spherical wind contributions have been omitted. 


In the case of the face-on narrow spiral, the spectral lines 
show no recognisable changes in the overall characteristics. 
However, other viewing angles show a very characteristic evolu¬ 
tion in the line shapes of the narrow spiral, with a progressively 
receding red or blue wing. This is explained as follows. For the 
higher CO transitions the emitting region, probing the most cen¬ 
tral regions of the spiral, will merely contain an incomplete por¬ 
tion of the first revolution of the spiral. Depending on whether 
the part of the spiral contibuting to the line flux is predominantly 
approaching or receding, emission is more concentrated in the 
blue or the red wing. It is likely that the most central regions 
of the binary system will show additional complexities, as this 
region harbours all the mechanics that form the spiral. These 
effects are not taken into account, and may cause additional ef¬ 
fects on the overall shape of the high CO transition lines. The 
shell spiral shows an evolution towards a triangular line shape. 
However, this is not dissimilar from the expected evolution of 
a homogeneous outflow, as the higher CO transitions probe the 
acceleration region of the wind. 


4.3. Spatial aspect 

In this section we present the spatial counterpart of the above re¬ 
sults. Generally, the two-fold nature of the wind is clearly recog¬ 
nisable in the PV diagrams. The spiral is manifested as a se¬ 
quence of emission bands, whilst the contribution of the spheri¬ 
cal wind is a large circular shape with a strongly enhanced cen¬ 
tral bar, representing the warm and dense regions near the star. 
This elongated feature approximately reaches up to the terminal 
velocity of the wind, and is visible in most PV diagrams below. 


4.3.1. Narrow spiral 

Inclination: The appearance of the narrow spiral is dependent 
on its orientation with respect to the observer, thus the veloc¬ 
ity channel maps will be sensitive to global angular transforma¬ 
tions. We recognise the central bar generated by the spherical 
wind in Fig. 0 The main feature, however, is the strongly peri¬ 
odic pattern seen in both PVl and PV2, generated by the spiral. 
The face-on view of the narrow spiral generates two identical 
PV diagrams, with a narrow column of periodic emission bands. 
The bands decrease in brightness for increased distance to the 
orbital plane (positioned at an angular deviation of 0 arcsec), 
as a result of the combined 1/r^ density dependence and local 
temperature. As the inclination increases, the primary contribu¬ 
tion to PV 1 gradually twists into an S-shaped feature, tightening 
around the central contribution, until it becomes a narrow, hori¬ 
zontal dumbell-like feature. In PV2 the evolution is different, be¬ 
cause of the strong asymmetry of the edge-on narrow spiral. Its 
evolution stretches the initial shape sideways, ultimately forming 
an ellipsiodal pattern of periodic emission bands for the highest 
inclination values. 

Mass contrast: As the mass contrast varies the emission 
contribution between the spiral and spherical winds changes ac¬ 
cordingly. In the case of E = 1/5, as seen in Fig.|9j the contribu¬ 
tion of the spherical wind is strongly enhanced, with a decreased 
relative emission of the spiral wind. However, the general struc¬ 
ture of the PV diagrams remains unchanged. 

Total mass-loss: In the high mass-loss regime the PV dia¬ 
grams show an overall increase in emission compared to the low 
mass-loss regime, as seen when comparing Fig. [TOl with Fig.O 
The apparent widening of the features of the spiral wind is due 
to non-negligible emission coming from the tails of the Gaus¬ 
sian density profile; emission which, in the low mass-loss case, 
barely contributed to the PV features. Optical depth effects keep 
the emission peak from increasing by the same factor as the wing 
emission. 
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Fig. 8: Position-velocity diagrams of a narrow spiral in a low mass-loss wind, in which the total amount of outstreaming mass is 
equally divided between both wind components, as seen from a range of different inclination angles. The top row shows PVl, the 
bottom row PV2. 
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Fig. 9: Position-velocity diagrams of a narrow spiral in the low 
mass-loss wind, in which the total amount of mass contained in 
the homogeneous outflow is five times the amount in the spiral, 
seen face-on and edge-on. The top row shows PVl, the bottom 
row PV2. 
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Fig. 10: Position-velocity diagrams of a narrow spiral in a high 
mass-loss wind for S = 1/5, seen face-on (left panels) and edge- 
on (right panels). The top row shows PVl, the bottom row PV2. 


4.3.2. Shell spiral 

The primary and spherical wind components possess identical 
spectral widths, and their general tendency to invariance under 
angular transformations translates into particularly similar PV1 
and PV2 diagrams. This strong charateristic is a reliable diag¬ 
nostic, assisting with the differentiation between the shell and 
narrow spiral. 

Inclination: Unlike for the spectral lines, the PV diagrams 
do reveal information on the inclination of the shell spiral, be it 
only for higher inclinations. Fig. [TT] shows that for low inclina¬ 
tions both PVl and PV2 are virtually indistinguishable. Higher 
inclinations show the appearance of gaps in the emission bands 
of PVl. The gaps appear at i=54 around -h and - 12 km/s and 
follow these bands as i increases, before reaching a verbal posi¬ 
tion at i=90. Under an iclination of i = 0°, any chosen PV slit 
(or axis along the plane of the sky) will be parallel to the orbital 
plane of the system. This means that the resulting PV diagram 
will be invariant to the particular choice of axis along which to 
construct the PV diagram. This explains why PVl and PV2 are 
identical. However, at an inclination of / = 0°, and when the axes 
are chosen to best show the asymmetry of the system, the orbital 
plane of the model will be perpendicular to one, whilst remain¬ 
ing parallel to the other axis. The axis parallel to the orbital plane 
thus again produces an identical PV diagram as before, whilst the 
other will enhance the features along the dimension perpendicu¬ 
lar to the orbital plane, which, in case of the shell spiral, contains 
the gaps. 

Mass contrast: The E = 1 / 5 case is displayed in Fig. [121 and 
shows how the relative contribution of the spherical wind inten¬ 
sifies as the contribution of the spiral wind fades. The general 
structure of the PV diagrams remains unchanged. 

Total mass-loss: A strong overall increase in emission, as 
seen when comparing Fig. [13] with Fig.[TT] The general structure 
of the PV diagrams remains unchanged. 

5. Discussion 

5.1. Constraining the geometry 

For our specific model parameters and our specific choice of PV 
axes, we find that PV2 readily provides the angular width and 
size of the gap between the emission bands, and thus the value 
of b. A slice through v = 0 kms“^ of PV2 allows to constrain 
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Fig. 11: Position-velocity diagrams of a shell spiral in a low mass-loss wind for E = 1, as seen from a range of different inclination 
angles. The top row shows PVl, the bottom row PV2 
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Fig. 12: Position-velocity diagrams of a shell spiral in a low 
mass-loss windforS = 1/5, seen face-on (left panels) and edge- 
on (right panels). The top row shows PVl, the bottom row PV2. 
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Fig. 13: Position-velocity diagrams of a shell spiral in a high 
mass-loss wind for E = 1, seen face-on (left panels) and edge-on 
(right panels). The top row shows PVl, the bottom row PV2. 
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the geometrical parameters in the orbital plane of the spiral. We 
refer to these diagrams as Slice Profiles (SPs) below. A number 
of such SPs are shown in Fig. [141 The top row exhibits the SP as 
they are, the bottom row shows the right half of the SP, overplot¬ 


ted with the geometrical model. We find as a general rule that the 
emission peaks systematically coincide with the density peaks. 
The angular distance between the spiral windings is thus readily 
available from the SPs of the PV diagrams. This distance can, if 
caused by bi nary interactions, dir ectly be related to local physics 
according to lKim & TaamI (1201 2l) via the relation 

/ 2 \ Invp 

Ara™ = J(yv.>+3Vpjx-^, (11) 

with (VV) the wind velocity, Vp the orbital velocity, and Tp the 
orbital radius of the primary. The orbital period Tp is given by 
the last term. 

The angular width of the emission bands, however, does not 
relate one-on-one to the spiral geometry and is thus more diffi¬ 
cult to deduce from actual observations. 

The plots labelled A in Fig.[T4]present the SP for a low mass- 
loss wind harbouring a narrow spiral seen edge-on, with E = 1. 
As expected, the spiral geometry tightly fits the emission charac¬ 
teristics. Additionally, a smooth bell-shaped background feature 
is visible, which is generated by the homogeneous component of 
the wind. Finally, the peak strength of the spiral-induced emis¬ 
sion spikes clearly changes as a function of offset. This is due 
to a combination of the model specific radial dependence of the 
spiral density and of the temperature law, as the latter defines the 
line-contribution regions. 

For more complex SPs, such as those for the shell spiral 
(labelled C), the emission peaks are not as narrow and smooth 
because the emission is smeared out over the offset dimension. 
Nevertheless, the presence of the spiral is still unmistakable. In 
this case, only the top portion of the peak, which ranges from ap¬ 
proximately the emission plateau strictly to the left of the spiral 
peak to the plateau strictly to its right, should be used to de¬ 
termine the geometry. Additionally, the bell-shaped background 
curve is still recognisable, but as a result of the nature of the shell 
spiral, the smeared-out emission somewhat conceals its shape. 

In the case of high mass-loss, the Gaussian tails of the den¬ 
sity distribution contribute considerably to the SPs, while optical 
depth effects prevent the highest density regions from augment¬ 
ing their emission by as much. This results in a decreased peak 
to background emission contrast. The emission from the tails 
broadens the emission zones significantly, which therefore no 
longer directly traces the width of the spiral arms. This effect 
is inherent to the mathematical properties of the Gaussian den¬ 
sity distribution. Sharp cut-offs of the spiral density would not 
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result in such a broadening. The bell-shaped emission feature 
of the smooth background of the A models has become much 
more prominent compared to the emission peaks generated by 
the spiral. Its shape has also changed to triangular. These two 
differences can be asigned to optical depth effects, which sup¬ 
press emission from high-density regions. 

It is necessary to note that these SPs should not be overinter¬ 
preted as their overall and absolute morphology is extremely sen¬ 
sitive to the specificities of the geometrical and radiative transfer 
model. Additionally, superimposed noise and diffraction effects 
will blur the small features in the SR 

5.1.1. Constraining the inclination 



Fig. 15: Left panel: A visual representation of the analytical 
model (green curve, Eq. [121) used to fit the emission features of 
the PV diagrams, applied to the reference model at an inclina¬ 
tion of 36°. The slope of the curve around v = Okms“^ (white 
dashed line) determines the characteristical skew of the diagram. 
Right panel: The relation between the measured skew of the PV 
diagram and the inclination angle. 


The dependence of the emission signature on the incination 
is characterised by the formation of an S-shaped feature in PV 1 
as inclination inceases. For a narrow spiral the S-shape is unmis¬ 
takable. The shell spiral exhibits this S-shape in the form of gaps 
in the emission bands. We constructed an analytical expression 
that fits the overall morphological trend of the inclination depen¬ 
dence of the PV 1 diagram. This relation is given by 


S(v) = -Ain 


B 



( 12 ) 


Exhibited as the wide green curve in Fig. o left panel, this ex¬ 
pression allows for the determination of the angle that the emis¬ 
sion feature makes with the velocity axis, around v = Okms“\ 
by using 


dS(v = 0) _4A 
Jv ” T’ 


(13) 


This is referred to as the skew of the diagram, and was deter¬ 
mined for a range of inclinations. We related the measured skew 
to the actual inclination angle to determine how they are related. 
The resulting curve (Fig. El right panel) shows a clear trend. 


5.2. Impact of the spiral structure on the total line strengths 

The signature of a spiral can quite easily be discerned in PV 
diagrams. However, the effect on the CO spectral lines is much 
less outspoken particularly so for extended (large a) spirals and 
high mass-loss rate systems. Here, we investigate the effect of a 
spiral structure on the frequency-integrated line strength of CO 
transitions. We discuss this integrated line strength relative to the 


case in which all of the gas is in a spherical outflow to assess the 
error we would make in constraining the mass-loss if we applied 
a model for a spherical wind to a system in which part (or all) of 
the material leaves the star in a spiraling outflow. 

Figure [T6| presents the total line flux normalised to the case 
of a spherical flow for six different CO transitions, from J=2-3 to 
38-37. The total flux ratio is plotted versus S. Models for a range 
of inclinations are presented at each value of E. The four basic 
configurations are labelled with different colours (blue for the 
shell spiral and red for the narrow spiral) and symbols (triangles 
for the low-M model and squares for the high-M model). The 
green horizontal bar represents a typical error margin of line- 
data calibrations, approximately 20 percent. 

For low values of E the bulk of the matter is in the spherical 
outflow and all ratios tend towards unity, as expected. The seem¬ 
ing failure of convergence to unity for high J levels originates in 
the fact that the elevated temperature of the spiral is not present 
in the spherical reference models. For the J=3-2 transition, opti¬ 
cal depth effects do not play a dominant role, safe for the high-M 
narrow spiral models. For these models the flux ratio drops be¬ 
cause of self-shielding. Note that, generally, the low-M models 
tend to have an increased flux ratio and the high-M models a 
decreased flux ratio. The former is an effect of density. If the 
material is concentrated in a spiral structure collisional excita¬ 
tions are more important causing the upper level of the transition 
to become overpopulated. This may increase the flux ratio by up 
to a factor of three. However, densities never reach values that 
strongly affect the populations of very high levels, with / > 30. 
Note also that the effect of inclination is not that important (com¬ 
pared to calibration uncertainties) for the low mass-loss case as 
the medium remains optically thin (in the line of sight). In the 
high-M case the line ratio drops, the more so if most of the ma¬ 
terial is in a narrow spiral structure (high E) and the inclination 
is high (approaching an edge-on view). In that case the flux may 
drop by up to a factor of ten. If the spiral fans out to well above 
the plane of symmetry (i.e. high a) the drop in flux is not that 
dramatic, simply because of less self-shielding compared to the 
narrow spiral. 

The models discussed here present a wide range in E, that 
may not all be realistic. The blue line in Fig. [16] represents the 
mass ratio E that was derived for the wind s t mctur e of the binary 
AGB system AFGL 3068 bv iKim & TaamI (l2Q12h . It shows that 
for this particular system one would not derive a significantly 
different total mass-loss rate if one would have applied a spheri¬ 
cal outflow model. 

We used the analytical expression derived by iDe Beck et al.l 
(|201Q|) , relating integrated line flux to mass-loss, to estimate 
the error on derived mass-loss rates provided prior knowledge 
on flux ratios. They showed that the analytically derived mass- 
loss scales as the integrated line flux to a certain exponent, 
M ~ [ J Fydvy, where y depends on the specific CO transition, 
and on the optical thickness of the spectral line. For the CO J=3- 
2 transition they found values for y to be 1.14 in the optically thin 
regime, and 0.61 in the optically thick regime. This translates, in 
the optically thin regime, into an uncertainty on the mass loss 
which is approximately a factor 1.3 greater than the uncertainty 
on the line flux (a flux uncertainty of a factor ten transflates into 
an uncertainty in deduced mass loss of a factor thirteen). For op¬ 
tically thick spectral lines, the uncertainty on the mass-loss is 
approximately a factor of 2.5 smaller than the uncertainty on the 
line flux (a flux uncertainty of a factor ten transflates to an uncer- 
tai nty in deduced mass lo ss of a factor four). Note that the results 
of jPe Beck et al.l (l201Ql) were derived by explicitly calculating 
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Fig. 14: A selection of SP for a S = 1 edge-on narrow spiral (A) in the low mass-loss regime, and (B) in the high mass-loss regime. 
(C) shows an SP with an identically parametrised shell spiral to (A). 
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Fig. 16: The ratio between the integrated flux of the spiral models, and the integrated flux of its ID equivalent, in which all material 
is located in a homogeneous outflow. In principle this flgure should exhibit identical amounts of red and blue markers, but a major 
overlap of both has resulted in an artiflcial bias favouring the red-labelled models. The green zone repr esents the instmmen tal 
calibration errors. The blue vertical line respresents the E value found in the hyrdodynamical simulations of iKim & TaamI (l2012h . 


the energy balance throughout the parameter study, as opposed 
to assuming a flxed temperature structure as in our models. 

It is worth noting that difficulties with the ID modelling of 
an abundance of spectral lines, like for example a combination 
of ground-based and unresolved PACS lines, are usually solved 
by forcing the model to have a peculiar and possibly somewhat 
misshapen temperature profile. The need for such actions might 
indicate the presence of a complex geometrical structure in the 
observed stellar wind. 


5.3. Single-dish simulations 

Spectral features of single-dish observations can be heavily in¬ 
fluenced by the telescope beam profile. To simulate how the 
intrinsic spectral lines of our models are affected by a beam 
profile, we covered our specific intensity images with a Gaus¬ 
sian Alter, centred onto the location of the central source. In 
Fig. [T7] the effect of Gaussian Alters with different full-width 
at half-maximum (FWHM) values on the overall shape of the 
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Fig. 17: The effect of single-dish telescope beam size on the spectral lines of three models is shown. The telescope beam size 
is simulated by a normalised Gaussian filter. The total flux of the spectral lines is therefore scaled by the normalisation factor 
1 /(cr V^), with FWHM = 2 ^/2ln{2cr). In effect, this means that the peak flux of the line (in the plot) scales with beam resolution. 
This is done to emphasise the evolution of the line shape. Left panel: shows the effect of a decreasing beam size on the face-on 
reference model, which creates a trident shaped profile. Middle panel: shows the effect on the edge-on reference model, which forms 
a very explicit double-peaked profile. Right panel: exhibits the beam-size effect on a low mass-loss, 2=1, face-on shell spiral. Here 
the double-peaked character of the spectral line is preserved. 


spectral lines is shown. The adopted values for the FWHM are 
^/20 for ^ € {1,2, 3,..., 10} of the total image size, which is 20 
arcsec. It is important to note that, for the clarity of the figures, 
the Gaussian filter was normalised. This operation scales the flux 
by the normalisation factor, allowing us to properly demonstrate 
the effect on the line shape without overlapping the lines. Plot A 
shows the effect of an ever narrower beam on a face-on narrow 
spiral. The prominent central feature remains, but as the beam 
size diminishes, a double-peaked profile forms around this cen¬ 
tral peak, created by the strongly resolved spherical wind. Plot 
B presents the effect of beam size on an edge-on narrow spiral. 
When fully resolved, the spectral line shows a very pronounced 
double-peak, which could be incorrectly identified as a detached 
shell. Finally, plot C shows the same effect for a face-on shell 
spiral. As a result of its similarity to a homogeneous outflow, the 
gradual evolution to a double-peaked profile is expected. How¬ 
ever, for an extremely resolved shell spiral, a small bump appears 
in the centre of the spectal line, which is some remaining numer¬ 
ical artefacts. 


5 . 4 . ALMA simulations 

The functions simobserve and the clean of CASA generated these 
synthetical ALMA obsevations. 


Simulation Parameters 

Pixel size 

0.04” 

Field size 

20” 

Peak flux 

Taken from LIME output 

Transition 

CO 3-2 (345.76599 GHz) 

Antenna configurations 

All available in cycle2 

Pointing 

Single 

Channel width 

1.15MHz, centred on rest freq. 

PWV 

0.913 mm 

Thermal noise 

standard 

Temperature 

269 K 

Integration time 

10 min on-source 


Table 3: The ALMA observation simulation specifications. 


In Fig. [TSl we present the effect of the seven different an¬ 
tenna configurations offered in ALMA Cycle2[3 on the face-on 
reference model, for which the fixed geometrical parameters can 
be found in Table[T] For these model-specific geometrical param¬ 
eters (combined with a fixed distance of 150 pc), we find that the 
resolution of the C34-1 configuration is too coarse. It only poorly 
samples the spiral structure. Additionally, the largest recoverable 
scale for the C34-1 configuration is large enough to nicely see 
the spherical wind. Combined, these effects show a strong bias 
towards the spherical wind. The C34-7 configuration has the op¬ 
posite effect, with a resolution high enough to detect the spiral 
structure, but a largest recoverable scale that is too small for a 
proper sampling of the spherical wind. Its flux contribution is 
completely lost. The C34-3 and C34-4 antenna configurations 
combine the advantages of both extremes, making them the op¬ 
timal setup with which to simulate observations of this spiral 
model. A direct consequence of using this best-fit configuration 
is that, as a general tendency, the overall structure of the PV 
diagrams remains relatively unaffected by the observation. We 
thus conclude that for this specific model, an extended config¬ 
uration with a maximum baseline of about 500m and a mini¬ 
mum baseline about 20m produces images of the higest quality. 
These baselines correspond to angular resolutions of approxi¬ 
mately 0.32-0.4 arcseconds, which, as expected, is similar to the 
angular size of the minimal length-scale of our models (the spi¬ 
ral arm thickness) of approximately 0.26 arcseconds. We would 
like to emphasise that these results are very model dependent, 
with the angular width of the smallest scale of the object being 
the most influential on the quality of the observations. The an¬ 
gular resolution A6 of a configuration with a maximum baseline 
L-max can be estimated with the following relation 


Lmax^ 

with y the observing frequency in GHz. Keeping the value of A6 
close to, but lower than the smallest expected angular size of the 
observed features of interest (in the case of the spiral structure 
the value (Xr/distance) should ensure a good quality observation. 

The compact array simulations are not shown. The resolution 
of the compact array is too low to be able to discern any useful 
features in the PV diagrams. 

^ http://almascience.eso.org/documents-and-tools/cycle-2/alma- 
technical-handbook 
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C34-1 C34-2 C34-3 c34-4 c34-5 c34-6 c34-7 



v[km/s] 

Fig. 18: A simulation of ALMA observations of the face-on reference model, with all the different extended antenna configurations 
offered during ALMA cycle two. The scaled colour bar used for the production of these images is identical to the one used to 
produce the intrinsic signatures, and is found in Fig. [3] 



A. intrinsic emission B. C34-1 C. C34-4 D. C34-7 

Fig. 19: The effect of three cycle 2 extended antenna configurations on the intrinsic emission of a number of face-on wind models, 
presented in panel A. Panel B shows the effect of the C34-1 configuration, with a minimum baslesine of 14.2m, a maximum baseline 
of 165.6m, and an angular resolution of 1.08 arcseconds. In panel C, the effect of the C34-4 configuration is shown. It is characterised 
by a minimum baseline of 20.6m, a maximum baseline of 558.2m, and a resolution of 0.32 arcseconds. Finally, shown in panel D 
is the influence of the C34-7 configuration, with a minimum baseline of 40.6m, a maximum baseline of 1507.9m, and a resolution 
of 0.12 arcseconds. 


Fig. [T9l presents the effect of the C34-l,4 and 7 antenna con¬ 
figurations on the PV diagrams of the face-on high and low mass- 
loss narrow spiral models (labelled by HN and LN respectively), 
and the high and low mass-loss shell spiral models (labelled by 
HS and LS respectively). An immediate conclusion from a pre¬ 
liminary inspection of the figures is that the lower resolution 
configuration provides considerably more information on the ob¬ 
ject than configurations with resolutions which are too high. The 
longest baseline configurations resolve out all the large-scale 
information, leaving only an exceptionally biased focus on the 
smallest scale, largest emission features. 

From these simulations we conclude that it is impossible, 
without any prior knowledge on the geometrical parameters of 
the spiral structure, to determine the optimal antenna configu¬ 
ration with which to observe the object. This being said, two 


scenarios can be thought of that can facilitate the configuration 
determination. Either hydrodynamical simulations predict a pre¬ 
ferred spiral arm thickness, inter-spiral-winding-distance, or re¬ 
lation between both, in which case an optimal (longest) interfer¬ 
ometry baseline can be calculated. Or theoretical models provide 
no such information, in which case it is advised to observe with a 
combination of configurations.To determine the best-fit antenna 
configuration the observed wind needs to be probed with at least 
the two most extreme maximum baselines. 


6. Conclusions 

Using the 3D radiative transfer code LIME, we have conducted a 
large-scale parameter study to investigate the effect of geometri¬ 
cal and global wind properties of AGB outflows containing spiral 
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structures that are embedded in a spherical outflow. Considering 
the CO v=0 J=3-2 transition, a singular spectral line generally 
conceals the dual nature of the wind. Only for limited combi¬ 
nations of parameters (inclination and spiral width) is its dual 
nature recognisable. However, when comparing different rota¬ 
tional transitions of CO, a characteristical evolution of the line 
shape is percieved as one progresses to higher transitions for the 
non-face-on narrow spiral models. The peculiarity of the (evolu¬ 
tion of the) resulting line shapes shows strong evidence for the 
presence of an embedded spiral. 

The LIME output also allowed us to investigate the 3D emis¬ 
sion of the CO v=0 J=3-2 transition throughout velocity space. 
We found that the best tools for analysing these images are the 
wide-slit PV diagrams, which intensify the correlated structural 
trends in the emission. Using these PV diagrams, we consistently 
found a strong signature of the spiral in the data, from which 
most of the geometrical properties can be recovered. It seems 
that one specific parameter, the mass- or density contrast, can¬ 
not be deduce by only analysing the emission, that is without 
resorting to detailed radiative transfer efforts. We compared the 
integrated flux of the 3D emission with an equivalent ID model 
(where all the material resides in the spherical outflow compo¬ 
nent) to explore the uncertainties such a ID misinterpretation 
brings about. We found that generally the errors on derived mass- 
loss do not exceed a factor of a few. Only in extreme cases can 
the difference between the embedded spiral flux and the exclu¬ 
sively spherical flux reach factors of up to ten, translating into 
uncertainties on derived mass losses exceeding an order of mag¬ 
nitude (for optically thin lines), or up to a factor of four (for 
optically thick lines). 

Finally, we have simulated cycle 2 ALMA observations of 
the model emission with CASA. In addition to the importance of 
recognising interferometric artefacts, we conclude that no pre¬ 
ferred antenna configuration exists, since the required resolu¬ 
tions depend strongly on the spiral geometry. To promote spiral 
detection the source of interest should be observed with a range 
of different maximum-baseline configurations. 
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